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Abstract 

Motivated by recent experiments on TmB4 (Phys. Rev. Lett. 101, 177201 (2008)), we examine 
the phase diagram of the Shastry-Sutherland antiferromagnet in an applied magnetic field in the 
limit of strong Ising anisotropy. In classical Ising limit, we demonstrate that the only fractional 
magnetization plateau is at 1/3 of the saturated magnetization. We study the perturbative influence 
of transverse quantum spin fluctuations, and present evidence that they can stabilize a narrow 1/2 
magnetization plateau. 



FIG. 1: The Shastry-Sutherland lattice. The exchange Ji acts between sites separated by the 
horizontal and vertical links, which the exchange J2 acts across the diagonal links. 



I. INTRODUCTION 

The Shastry-Sutherland antiferromagnetP is a rare example of an exactly soluble frus- 
trated spin model in two spatial dimensions. Its ground state is a direct product of near- 
neighbor spin-singlet pairs, and there is an energy gap to all excitations. Much theoretical 
and attention has been lavished on this model after the discoverj^^ of its experimental 
realization in SrCu2(B03)2. 

This paper will examine the phases of the Shastry-Sutherland antiferromagnet for the case 
of exchange interactions with a strong Ising anisotropy, and in the presence of an applied 
magnetic field. This is motivated by a recent experimental studj^ of TmB4, in which the 
magnetic Tm ions are believed to be Ising-like. 

The experimental studies on the compound TmB4 determined both on its phase diagram 
as a function of temperature, and the magnetization diagranP as a function of the external 
magnet field. In the magnetization diagram, as the external magnetic field increased from 
0, magnetization plateaus with M/Msat = 1/7, 1/8, 1/9... are seen, followed by a major 1/2 
magnetization plateau before reaching saturatiorP. 

In this paper a classical Ising model was first used in Monte Carlo calculation to simulate 
the magnetization process. Based on the simulation result, a more detailed analysis was 
done on the ground state phase diagram for Ising-like Shastry-Sutherland (SS) lattice. Going 
beyond the classical Ising limit, we accounted for the quantum effects of the transverse spin 
interaction by the effective Hamiltonian method of Sen, Wang and Damle^; this is valid in 
the limit of large spin S. We will show that this effective Hamiltonian does contain terms 
which stabilize a 1/2 magnetization plateau. 

The theoretical models in this paper were all based on the two-dimensional lattice, with 
magnetic ion located on each site of the SS lattice. The magnetic ion interacts only with the 
nearest neighbors through exchange interactions, with square bond strength Ji and diagonal 
bond J2 (see Fig. [T]). 

While this work was being completed, we learnt of the work of Meng and WesseP; where 
their results overlap with ours, they are in agreement. 
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FIG. 2: The magnetization process when J2 take the three different values: J2 
4Ji. 



Jl, J2 — 2 Jl, J2 



II. THE SSL GROUND STATE PHASE DIAGRAM IN THE ISING MODEL 

A naive approach to the problem is the classic spin-1/2 Ising model. Spins, which can 
take two different values (up or down), are located at each corner of the SS lattice, with 
anti-ferromagnetic exchange interaction Ji > and J2 > between the nearest neighbors. 
To simulate the magnetization process, the Monte Carlo Method was used in the simulation. 
The simple Hamiltonian was taken as: 

H = Jl aiaj + J2 <Ji<Jj — fiB a, (2.1) 

Jl(i,j) J2(i,j) i 

where the first summation runs over all the square bonds and the second runs over all the 
diagonal bonds, a denotes the spins with cTj = ±1, and fiB term denotes the energy gained 
from the external magnetic field. 

The detailed Monte Carlo process was done in the following way. The lattice size was 
chosen to be 18 x 18, which is a common multiple of both 2 and 3 (the reasons will follow), 
with periodic boundary condition at the four sides. Periodic boundary condition in SS lattice 
requires the length to be a multiple of 2. In the Metropolis algorithm, iteration times before 
reaching equilibrium and after reaching equilibrium were all set to be 100,000 times, the 
temperature was set very close to zero, and the bond strength Ji was set to one. 

The results of the simulation are shown in Fig. [2| where different lines represent different 
values of J2/J1 ratio. As can be seen from the figure, for all values of the J2/J1 ratio, there's 
always a very large 1/3 magnetization plateau, and no 1/2 plateau is observed. It will be 
shown later in this section that no 1/2 plateau is possible at zero temperature in the Ising 
model. Similar simulation results were also be obtained when lattice size was chosen to be 
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FIG. 3: The spin configuration pattern obtained in the one-third magnetization plateau. 



12 or 24 (again, both as integer multiple of 2 and 3). It is important to note that in this 
simulation, the choice of the lattice size matters, and a wrong choice of lattice size could 
lead to different conclusions. To understand the choice of lattice size, the resulting spin 
configuration in our 1/3 magnetization plateau is shown here in Fig. |3} 

In Fig. [3} empty dots represent spin ups while black dots represent spin downs. In the 
figure, 1/3 of the total spins are pointing down while others pointing up. Also, every spin 
down is surrounded by spin ups both through Ji bond and J2 bond. This is the key feature 
of this spin configuration, which makes its energy the lowest among all that have 1/3 of 
spins pointing down (c.f. Eq. ( 2.1[ )). Due to the period of 3 in this spin configuration (a 
repetition of three rows in the figure), the lattice size has to be a multiple of 3 to make 
this configuration possible. Note that the degeneracy of this spin configuration has a very 
limited value of 6, compared to the large degeneracy of other types. From the figure, we can 
now write down the Hamiltonian for this particular spin configuration (Hamiltonian for the 



1/3-plateau phase) by using Eq. (2.1): 



H = -\NJ^-\NJ2-\N^^B 
3 3 



(2.2) 



where represents the total number of lattice sites, and other parameters are the same as 



Eq. (2.1) 



We are now at the position of giving out a full ground state phase diagram. The proposed 
ground state phase diagram is shown here in Fig. |4j Different phases shown are the different 



zero-temperature ground states of the general SS lattice with Hamiltonian of Eq. (2.1). As 



shown in the figure, when the external field is zero, we have the Neel State (opposite spins 
on square bonds) when J2/J1 ^ 1, and the Ising Dimer State (opposite spins on diagonal 
bonds) when J2/J1 ^ 1. As the external field increases, the 1/3 plateau phase will arise, 
giving us a positive total magnetization. When the external field is large enough, the All 
Spin- up State will show up, giving us the magnetization saturation. 

Further reflection explains the existence of 1/3 plateau: if we rewrite the Hamiltonian 
as a function of the ratio between spin-down-number and total-spin-number (denoted by 
p), the 1/3 plateau phase will have the largest p value among all states that have the key 
feature (all the spin downs are surrounded by spin ups through both Ji and J2 bonds). That 
is to say, the 1/3 plateau state will have the largest number of spin downs while keeping 
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FIG. 4: The ground state phase diagram for SS lattice, with J2/J1 and /iB/Ji being two inde- 
pendent variables. The boundaries are obtained simply by comparing the Hamiltonian of different 
phases using Eq. (2.1). 



the key feature. This will make the 1/3 plateau phase an extreme. Different extremes of 
the Hamiltonian are obtained at p = 0(all spin-up state), p = 1/3(1/3 plateau phase) and 
p = 1/2 (Neel state and Ising dimer state), which are exactly the different phases shown in 
Fig.i 

To demonstrate the reliability of the proposed ground state phase diagram, an exhaustive 
computational simulation was done on a 4 x 6 lattice, by going through all the possible spin 
configurations and picking out the ones with the lowest energy (altogether 2^^ ^ 1.68 x 10'' 
configurations). The result of the simulation agreed completely with Fig. |4| 

It's also very instructive to look back and compare Fig. |4]with the magnetization result 
obtained in Fig. [2j The two figures agree well at the critical points when phase transition 
occurs from the 1/3 plateau to the all spin-up state, with J2/J1 value being 1, 2, 4, and 
pB/Ji value being 5, 6, 8, respectively. 

To conclude, as a direct application of the obtained phase diagram, we'll present here a 
short proof showing that no 1/2 plateau ground state is possible at zero temperature in the 
Ising model. This is achieved by showing that there's always a lower energy state than any 



possible states with a 1/2-saturation magnetization. First, in view of Eq. (2.1), it's not hard 
to get the lowest energy expression for states with a 1/2 magnetization: 

H = -^NpB (2.3) 

In the above equation, by minimizing the energy expression, both Ji and J2 terms ac- 
cidentally vanished, leaving only the pB term. The minimizing process could be easily 
understood as follows: for states with a 1/2-saturation magnetization, the total number of 
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spin downs would be 1/4N{N being the total site number), so that the maximum number 
of antiparallel square-bond pairs would be 1/4N x 4 = N. Noting that there're altogether 
2N square-bonds, the antiparallel pairs and parallel pairs would cancel each other in the Ji 
term. Similar analysis applies to the J2 term. 

The energy expression for the All Spin-up State is also needed for the proof: 



H = 2NJi + y J2 - NfiB 



(2.4) 



Now, by direct comparison of the three energy expressions (Eqs. 2.2, 2.3| and 2.4), we see 



that either the All Spin-up State (Eq. ( |2.4 )) or the 1/3 Plateau Phase (Eq. (2.2)) have lower 
energy than the minimum energy we obtained for states with 1/2-saturation magnetization 



(Eq. (2.3)) throughout the phase diagram. That is to say, states with 1/2-saturation magne- 
tization could never have the lowest energy. Thus no 1/2 magnetization plateau is possible 
in the Ising limit at zero temperature. 



III. PERTURBATIVE INFLUENCE: TRANSVERSE QUANTUM FLUCTUA- 
TION 



As the classical Ising model does not give us an understanding of the major 1/2 plateau 
observed in the experiment,'^ improvements were made to our Ising model Hamiltonian. In 
this section, we focus on the perturbative effects of the transverse fluctuation on the SSL 
spin multiplets with large easy axis anisotropy: 

H = - DY^iS!)' + J, S!S^ + J2 sts^-fiBYs: 



Ji Y iS!Sj + sfs]) + J, Y (^r^J + ^'^J) 



(3.1) 



In the above equation, D is the strength of the easy axis anisotropy along the z axis. Si are 



the spin multiplets, Ji and J2 are the bond strengths same as in Eq. (2.1), and the p term 



was introduced to take care of the transverse interaction, with parameter p adjustable to 
count for the asymmetry between z axis and the x,y axis. 

As argued in Ref. |5l since the problem turns out to be a strong Ising anisotropy problem 
(see Ref. a J/D expansion of the Hamiltonian will give us the leading quantum effects. 
For the perturbation process, following the method introduced by Ref. |5l we first split the 
Hamiltonian as H = Hq + H' , in which 



^0 



-DYis:f + Ji Y sts] + j2 Y sts]-iiBYs! 



(3.2) 



J2{i,j) 
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(3.3) 



J2{i,j) 



Note that H' is actually the p term in Eq. (3.1). When D is large, the ground state 



manifold of Hq is just the Ising limit we adopted in the previous section, where the spins 
are all maximally polarized along the z axis (so that = Sdi). This ground state manifold 
is referred to as the Ising subspace. Degenerate perturbation theory will now be applied by 
treating H' as a perturbation on this Hq ground state manifold (or, the Ising subspace). 

It's not hard to check that both the first-order perturbation terms and the second-order 
off-diagonal perturbation terms vanish, leaving only the second-order diagonal terms: 



\{0\H' 



k^G 



(3.4) 



where G denotes the Ising subspace, denotes any state in the Ising subspace, k denotes 
any state outside the Ising subspace, and Eq and Ek are the energies for the corresponding 
states. For further convenience, rewrite the above expression in the following way: 



HP) 



E 



\m'\ki] 

Eq — Ek 



E 



En — E, 



k2 



(3.5) 



where ki are the states that, compared with states in the Ising subspace, have different spins 
only at the two ends of one single square bond, while ^2 having different spins only at the 
two ends of one single diagonal bond (it is not hard to see that ki and k2 will include all the 



states with non-vanishing value in Eq. (3.4)). 

After a somewhat lengthy calculation of the above equation by calculating the two terms 
separately (note that in the calculation, the two terms would have different energy denomi- 
nators), and by performing a small- J/D expansion for both terms, we finally arrive at the 
resulting perturbation term: 



1 - aiCTj 



+ 



2D{2S-1) 2 

J2S{(yiOid + g-jg-jrf) 
2D(25- 1) 



1 + 



■h 



2D(2^-1) ' 2D(2^- 1) 



+ 



2D(2^ 



J2{id) 



2D 2D{2S-1) 



(3.6) 



and the effective Hamiltonian in the Ising subspace was obtained as: 

H = Hq + H^^^ 



(3.7) 



In Eq. (3.6), S is the absolute value of the spin multiplet, CTj = ±1 are the spin-ups or 
spin-downs as in the previous Ising limit approach, whereas Hi = ctijO'j is the local field, 
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FIG. 5: The possible "new" phase diagram with the new effective Hamiltonian. In the diagram, all 
the boundaries between different phases are obtained using the effective Hamiltonian. Note that 
in the figure, there's now a split between the "All Spin-up State" and the "1/3 Plateau Phase", 
where a state with magnetization being 1/2 of the saturation value could be proved to have lower 
energy than both of the two states, and thus provide a region for a possible "1/2 Plateau Phase". 



with = 1 for aj's square-bond neighbors and aij = otherwise; and aid 
aij = 1 for (Ji's diagonal bond neighbor and aij = otherwise. 



dijCTj, with 



A. The Nev^ Phase Diagram 

With the new effective Hamiltonian obtained above, we can now propose a new ground 
state phase diagram for the system. Note that this is far from a rigorous proof of the diagram, 
as the Hamiltonian is much more complex than in the Ising model, and there's now no direct 
way to get the global energy minimum. Also, the exhaustive simulation is impossible in this 
case. Yet as the new Hamiltonian is based on the original Ising Hamiltonian, it's reasonable 
to start from the ground states in the original phase diagram. 

In Fig. |5| we directly adopted the previous four ground states without rigorous proof. 
The boundaries in the figure are all obtained using the new Hamiltonian, where we've set 
the magnitude of the spin multiplet S' to be 6 (as in the experimeniP) , the Ji/D value to be 
1/5, and the parameter p to be 1. For the convenience of direct comparison with the original 
phase diagram, the unit of the vertical axis was divided by a factor of 5*. Comparison with 
the original phase diagram (Fig. |4]) shows their similarities, which should be the case. 

Of course, minor corrections are necessary for the new phase diagram, especially at the 
boundaries between the four "old" ground states; here presented in Fig. [5] is only a correction 
we found at the boundary between "All Spin-up State" and the "1/3 Plateau Phase". As 
could be seen in the figure, the boundary now split up, and a new state with magnetization 
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FIG. 6: The spin configuration pattern of the "new possible ground state" found. As could be seen 
in the figure, 1/4 of the total spins are pointing down, making the total magnetization 1/2 of the 
saturation value. 
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FIG. 7: The phase diagram obtained with a more dramatic p value. The parameter -p is changed 
to 5 in the figure, to make the possible "1/2 Plateau Phase" more obvious to see. Note that in the 
figure, the boundaries began to cross each other, which means that the boundary behavior of the 
phase diagram is more complex than before, and further analysis are needed. 



being 1/2-saturation is found to have lower energy. Its detailed spin configuration pattern 
is showed here in Fig. |6] (again, black dots represent spin downs while empty dots ups). 
Although it's impossible to prove that this configuration has the lowest energy, this inter- 
esting finding at least gives a region where the "1/2 magnetization plateau" is possible. If 
we redraw a diagram of magnetization against the external field, before the magnetization 
reach its saturation value, a 1/2 plateau may occur. To make this region more dramatic and 
easy to see in the phase diagram, the parameter p was tuned to be 5 in Fig. [7j 

At the end, it's also reasonable to further suspect that the small fractional magnetization 
plateaus with MjM^at = 1/7,1/8,1/9... found in the experimentpl could also be possibly 
explained by similar arguments. As the boundary between "Neel State" and "1/3 Plateau 
Phase" start to split up, fractional plateau phase may have lower energy and arise from the 
effective Hamiltonian correction. 
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IV. CONCLUSION 



In this paper, a thorough analysis was done for the two dimensional SS lattice with 
the Ising model, giving us the full ground state phase diagram (see Fig. |4]). The phase 
diagram is consistent with both the Monte Carlo simulations and the exhaustive simulation 
for a small system. As shown in the diagram, except for the Neel State, Ising Dimer State 
and the All Spin-up State, only a 1/3 magnetization plateau phase could be found, and 
no 1/2 magnetization plateau is possible in the Ising limit. Further study was done by 
considering the perturbative effects of the transverse fluctuations, followed by J/D strong 
anisotropy expansion. The new effective Hamiltonian showed minor corrections to the phase 
diagram, and more importantly, it gave a possible region in the phase diagram where 1/2 
magnetization plateau may occur. Further analysis could be done in the future to confirm 
this "1/2 plateau region", or to probably explore corrections at other boundaries, in the 
hope of finding other small "fractional plateau phases" in this model. We hope that the 
work done would be useful for further understanding of the behaviors of Tmi?4, and for the 
further study of all the SSL quantum magnets. 



Acknowledgments 



This research was supported by the NSF under grant DMR-0757145. 



B. S. Shastry and B. Sutherland, Physica B 108, 1069 (1981). 

H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, K. Onizuka, M. Kato, K. Kosuge, C. 
P. Slichter, T. Goto and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 

H. Kageyama, M. Nishi, N. Aso, K. Onizuka, T. Yosihama, K. Nukui, K. Kodama, K. Kakurai, 
and Y. Ueda, Phys. Rev. Lett. 84, 5876 (2000). 

K. Siemensmeyer, E. Wulf, H.-J. Mikeska, K. Flachbart, S. Gabani, S. Mates, P. Priputen, 
A. Evdokimova, and N. Shitsevalova, Phys. Rev. Lett. 101, 177201 (2008). 
A. Sen, F. Wang, and K. Damle, |arXiv:0805.2658| (2008). 
Z. Y. Meng and S. Wessel, |arXiv:0808.3104 (2008). 



10 



